PERSISTENCE AND PERMANENCE OF MASS- ACTION AND 
POWER-LAW DYNAMICAL SYSTEMS 
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Abstract. Persistence and permanence are properties of dynamical systems that describe the 
long-term behavior of the solutions, and in particular specify whether positive solutions approach the 
boundary of the positive orthant. Mass-action systems (or more generally power-law systems) are 
very common in chemistry, biology, and engineering, and are often used to describe the dynamics in 
interaction networks. We prove that two-species mass-action systems derived from weakly reversible 
networks are both persistent and permanent, for any values of the reaction rate parameters. More- 
over, we prove that a larger class of networks, called endotactic networks, also give rise to permanent 
systems, even if we allow the reaction rate parameters to vary in time. These results also apply to 
power-law systems and other nonlinear dynamical systems. In addition, ideas behind these results 
allow us to prove the Global Attractor Conjecture for three-species systems. 
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1. Introduction. Determining qualitative properties of solutions of dynamical 
systems arising from nonlinear interactions is generally a daunting task. However, a 
relevant mathematical theory that pertains to biochemical interactions obeying mass- 
action kinetics has been developed over the last 40 years, starting with the seminal 
work of Fritz Horn, Roy Jackson and Martin Feinberg, [HI |TTJ [T2] . Generally termed 
"Chemical Reaction Network Theory" , this theory establishes qualitative results that 
describe the surprisingly stable dynamic behavior of large classes of mass-action sys- 
tems, independently of the values of the reaction rate parameters [8]. This fact is 
especially useful since the exact values of the system parameters are usually unknown. 

Here we focus on the properties of persistence and permanence for mass-action 
systems, and the more general power-law systems. A dynamical system on K™ is 
called persistent if no trajectory that starts in the positive orthant has an w-limit 
point on the boundary of R" > an d is called permanent if all trajectories that start 
in the positive orthant eventually enter a compact subset of R" . Persistence and 
permanence are important in understanding properties of biochemical networks (e.g., 
will each chemical species be available indefinitely in the future), and also in ecology 
(e.g., will a species become extinct in an ecosystem), and in the dynamics of infectious 
diseases (e.g., will an infection die off, or will it infect the whole population). 

In the context of biochemical networks we formulate the following conjecture: 

Persistence Conjecture. Any weakly reversible mass-action system is persistent. 

Here a weakly reversible mass-action system is one for which its directed reaction 
graph has the property that each of its connected components is strongly connected^ 
[8]. A version of this conjecture was first mentioned in [9^ 

In this paper we prove the Persistence Conjecture for two-species networks. More- 
over, we introduce a new class of networks called endotactic networks, which contains 
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the class of weakly reversible networks, and we show that endotactic two-species in- 



variable mass action systems are permanent (Theorem 6.4 1. 

For two-species systems, this is a stronger statement than the Persistence Con- 
jecture, in several ways. First, permanence here means that there is a compact set 
K in the interior of the positive orthant such that all forward trajectories of positive 
initial condition end up in K. Therefore, our result implies both strong persistence 
and uniform boundedness of trajectories for two-species K-variable mass-action sys- 
tems. Second, our mass-action system is n-variable, meaning that the reaction rate 
constants are allowed to vary within a compact subset of (0,oo). This makes our 
result applicable to a larger class of kinetics than just mass-action. Finally, the class 
of endotactic networks, explained in detail in section [4] strictly contains the class of 
weakly reversible networks. We conjecture that endotactic K-variable mass action 
systems are permanent for any number of species. 

On the other hand, note that the Persistence Conjecture for three or more species 
remains open. 

The Persistence Conjecture is also strongly related to the Global Attractor Con- 
jecture, which was first stated over 35 years ago [TT], and is one of the main open ques- 
tions in Chemical Reaction Network Theory. It is concerned with the global asymp- 
totic stability of positive equilibria for the class of "complex-balanced" [TH [TTJ [5J [5] 
mass-action systems. It is known that such systems admit a unique positive equi- 
librium cs within each stoichiometric compatibility clasfj^j S. Moreover, each such 
equilibrium admits a strict Lyapunov function and therefore is locally asymptoti- 
cally stable with respect to £ [SI HI]. However, the existence of this Lyapunov function 
does not guarantee that cs is a global attractor, which is the object of the following 
conjecture: 

Global Attractor Conjecture. Given a complex-balanced mass-action system and 
any of its stoichiometric compatibility classes S, the positive equilibrium point cs is 
a global attractor on mi(£). 

Trajectories of complex-balanced mass-action systems must converge to the set 
of equilibria [15] . This makes the Global Attractor Conjecture equivalent to showing 
that any complex-balanced system is persistent, i.e. no trajectory with positive initial 
condition gets arbitrarily close to the boundary of the positive orthant. Therefore, 
since a complex-balanced mass-action system is necessarily weakly reversible [5J, a 
proof of the Persistence Conjecture would imply the Global Attractor Conjecture. 
In general, both conjectures are still open; the Persistence Conjecture has not been 
proved previously in any dimension, but recent developments have been achieved to- 
wards a proof of the Global Attractor Conjecture. For instance, Anderson [1] and 
Craciun, Dickenstein, Shiu and Sturmfels [5] showed that vertices of the stoichio- 
metric compatibility class S cannot be w-limit points and that the Global Attractor 
Conjecture is true for a class of systems with two-dimensional stoichiometric compat- 
ibility classes. Anderson and Shiu proved that for a weakly reversible mass-action 
system, the trajectories that originate in the interior of £ are "repelled" away from 
the codimension-one faces of S. This allowed them to prove that the Global Attractor 
Conjecture is true if the stoichiometric compatibility class is two-dimensional [2J. 

2. Definitions and notation. A chemical reaction network is usually given by 
a finite list of reactions that involve a finite set of chemical species. For example, a 



reaction network involving two species A\ and A2 is given in Figure 2.1 



3 A stoichiometric compatibility class is a minimal linear invariant subset; see Definition 2.2 
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Fig. 2.1. A two-species reaction network. 



We define the functions ca 1 (i) and ca 2 (t) to be the molar concentrations of species 
A\,Ai at time t. The changes in the concentrations are dictated by the chemical 
reactions in the network; for instance, whenever the reaction Ai + A2 — > 2Ai occurs, 
the net gain is a molecule of Ai, whereas one molecule of A2 is lost. 

2.1. Chemical reaction networks. Here we recall some standard terminology 
of Chemical Reaction Network Theory (see [12] ). In what follows the set of non- 
negative, respectively strictly positive real numbers are denoted by K> and K>o . 
For any integer n > 1 we call K™ the positive orthant and M™ the closed positive 
orthant. For an arbitrary finite set / we denote by Z> and M> the set of all formal 

sums a = otii where are nonnegative integers, respectively reals. The support 
iei 

of an element a € M. 1 is supp(a) = {i 6 / : a.% 7^ 0}. 

Definition 2.1. ^4 chemical reaction network is a triple (S,C,TZ), where S is 
the set of chemical species, C C Z> is the set of complexes (i.e., the objects on left 
or right side of the reaction arrows), and 1Z is a relation on C, denoted P — > P' , 
representing the set of reactions in the network. Here P is called the source complex 
and P' is called the target complex of the reaction P — > P' . Moreover, the set TZ must 
satisfy the following three conditions: it cannot contain elements of the form P — > P; 
for any P 6 C there exists some P' G C such that either P — > P' or P' — > P; and the 
union of the supports of all P G C is S . 

The second condition in Definition |2.1| guarantees that each complex appears in at 
least one reaction, and the third condition assures that each species appears in at 
least one complex. For the reaction network in Figure |2.1[ the set of species is S = 
{Ai, A2}, the set of complexes is C = {2Ai, A\ + A2, 2A2} and the set of reactions is 
TZ = {2Ai ^ Ai + A 2 ,Ai + A 2 ^ 2A 2l 2A 2 ^ 2Ai], which consists of 6 reactions, 
represented as three reversible reactions. For convenience, we refer to a chemical 
reaction network by specifying TZ only, since it encompasses all the information about 
the network. 

By a slight abuse of notation we may view all complexes as (column) vectors of 
dimension equal to the number of elements of S, via an identifica tion given by a fixed 



ordering of the species. To exemplify, the complexes in Figure 



Ai + A 2 = 



and 2A 2 = 



2.1 



are 2Ai = 



For any reaction P — > P' we may thus consider the 

vector P' — P, called the reaction vector of P —> P' . 

Definition 2.2. The stoichiometric subspace S of TZ is defined as 



S = span{P' - P I P -> P' € TZ}. 

Suppose TZ has d species, fix an order among them and c(t) eR s = R d denote the 
vector of molar species concentrations at time t. We refer to c(t) as the concentration 
vector of TZ at t. As we will see soon, for all t > 0, the concentration vector c(t) 
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is constrained to the stoichiometric compatibility class of c(0), a special polytope in 
K> , denned below. 

Definition 2.3. Let c° G R d . The polytope (c° + S) ("I M> is called the stoichio- 
metric compatibility class ofc°. 

The reaction network 7Z can be viewed as a directed graph whose vertices are 
given by complexes and whose edges correspond to reactions of 1Z. The connected 
components of this graph are called linkage classes of 1Z. 

Definition 2.4. A reaction network 1Z is called reversible if P' — » P G 1Z 
whenever P — > P' 6 TZ. 

Definition 2.5. A reaction networklZ is called weakly reversible if its associated 
directed graph has strongly connected components. 

2.2. K-variable mass-action systems. Throughout this paper, the differential 
equations that govern the evolution of the concentration vector c(t) will be given by 
K-variable mass-action kinetics. The next definition clarifies this notion; to state it, 
we need the following notation: given two vectors u — ^""^ u s s and v = Vju s s in 

ses sS5 
M> , we denote u v — J|(u s )"', with the convention 0° = 1. 

sGS 

Definition 2.6. A n-variable mass-action system is a quadruple (S,C,7Z, k) 
where (S,C,1Z) is a reaction network and k : R>o — > {rj, l/^) 7 ^ for some r\ < 1 is a 
piecewise differ entiable function called the rate constants function. Given the initial 
condition c(0) = c° G M^ , the concentration vector c(t) is the solution of the K- 
variable mass-action ODE system 



;(*)= Yl ^P>(t)c(t) p (P'-P). (2.1) 



p^P' 



The /{-variable mass-action kinetics is a natural generalization of the widely 
used mass-action kinetics, where the rate constant functions «(•) are simply positive 
constants. Therefore any results we prove for K-variable mass-action systems are also 
true for mass-action systems. 



Note that integration of (2.1) yields 



c(t) = c° + £ QT K P ^ P ,(s)c(s) P d S ^ (P' - P), 

and therefore c(t) G c(Q) + S for all times t > 0. Moreover, it is easy to check that 
c{t) G R>o; indeed, all the negative terms in the expression of Cj(t) coming from the 



right hand side of (2.1 ) contain Cj(t) as a factor and vanish on the face c, = of 



We conclude that (c u + 5) n E> is forward invariant for (2.1). Similarly one can 
notice that (c° + S) PI M> ^ s a l so forward invariant for (2.1 ). 



2.3. Persistence and permanence of dynamical systems. We review some 
more vocabulary ( |101 [TBI, H^] ). 

Definition 2.7. A d-dimensional dynamical system is called persistent if for 
any forward trajectory T(cq) = {c(t) = (xi(t), . . . ,Xd(t)) \ t > 0} with positive initial 
condition c G M> we have 

lim inf Xi (t) > for all i G {1, . . . , d}. 

t—*oo 



5 



Note that some authors call a dynamical system that satisfies the condition in 
Definition 2.7 strongly persistent [19] . In their work, persistence requires only that 
lim sup^go Xi (t) > for all i e {1, ...,d}. 

Going back to Definition |2.7| persistence means that no forward trajectory with 
positive initial condition approaches the coordinate axes arbitrarily close. Note that a 
dynamical system with bounded trajectories is persistent if the open positive quadrant 
is forward invariant and there are no w-limit points on <9R> ; recall the definition of 
w-limit points below. 

Definition 2.8. Let T(c ) = {c(t) \ t > 0} denote a forward trajectory of a 
dynamical system with initial condition cq 6 Kl>o- The w-limit set of T(cq) is 



lim^ T(co) = {I £ M. d | lim n _ i . 00 c(t n ) — I for some sequence t n — > oo}. 

The elements o/lim w T(co) are called u -limit points o/T(co). 

Definition 2.9. A d-dimensional dynamical system is called permanent on a 
forward invariant set D C R> if there exists e > such that for any forward trajectory 
T(cq) with positive initial condition cq € D we have 



e < liminf Xi(t) and lim sup X;^) < 1/e for all i £ {1, 

t-^oo t-yoo 



In other words, a dynamical system is permanent on a forward invariant set 



D C 



if there is a compact set K C 



such that any forward trajectory with 



positive initial condition cq € D ends up in K. 

In the case of a K-variable mass-action system, where the trajectories are confined 
to subspaces of M. d , it is meaningful to require that the system is permanent on each 
stoichiometric compatibility class, as we do in the following definition. 

Definition 2.10. A n-variable mass-action system is called permanent if, for 
any stoichiometric compatibility class S, the dynamical system is permanent on E. 
Clearly, a permanent k- variable mass-action system is persistent. 

3. An illustrative example. Before introducing endotactic networks and prov- 



ing their persistence (Theorem 6.2 1, we take an intermediary step and discuss an ex- 



ample which illustrates some key points in our strategy for the proof of Theorem 6.2 
Consider the following K-variable mass-action system: 



k\ k-2 &3 

2X X X ^2X + Y 

fc-l fc-f! 



(3.1) 



where fcj = fcj(i) G (77,1/77) for all i € {— 3, — 2, . . . , 3} and all t > 0. Let c(t) = 
(x(t),y(t)), where x{t) and y(t) denote the concentrations of species X and Y at time 
t > 0. Then (x(t),y(t)) satisfy the system of differential equations (3.2): 



x(t) 

m 



-[k 2 x(t)-k- 2 y{t)} 



-[k 3 x(t)-k_ 3 x(t) 2 y{t)} 



(3.2) 

Let us fix an initial condition c(0) for (3.2). We want to construct a convex polygon V 



such that conv(V) is forward invariant for the dynamics given by (3.2) and contains 
c(0), where conv{V) is the union of V and its interioiQ According to a theorem of 
Nagumo (see [14j 0]), this is equivalent to requiring that the flow satisfy the sub- 
tangentiality condition, which in our case simply means that the vector c(to) points 



4 In this case we will simply say that V is forward invariant for (3.2 I. 



G 




towards the interior of V whenever c(io) S "P- In example (3.1), the flow has three 
components along reaction vectors corresponding to the three pairs of reversible reac- 
tions (see (3.2)). We would like to construct V such that each component of the flow 
satisfies the sub-tangentiality condition; this will clearly imply that the aggregate flow 
(3.2) also satisfies the sub-tangentiality condition and therefore V is forward invariant 
for |3.2). 

Following up on the preceeding observation, we next try to characterize the poly- 
gons that are forward invariant for one component of the flow, i.e., corresponding to a 
single reversible reaction. Let us fix t for a moment. The first reaction of (3.1 ) drives 

_2] 

with rate k\{t)x{t) 2 — k—i(t)y(t). The sign of the rate 



the flow along the vector 



1 



gives the direction of the flow, 



-2" 




' 2 " 


1 


or 


-1 



according to which side of the curve 



y = (ki(t)/k_i(t))x 2 the phase point c(t) belongs to. The coefficient ki(t)/k-i(t) 
is unknown, but we know it is bounded, fci(i)/fc_i(i) £ (n 2 ,l/ri 2 ). Therefore, (see 
Figure |3~lj (a)) the component of c(t) in ( |3.2[ ) corresponding to the first reaction is 

if c(t) G {(x,y) € M> | y > {1/S)x 2 } and in direction 



in direction 



-1 



1 



if 
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c(t) 6 {(x,y) € M>o I 2/ < <5^c 2 } , where <5 = rj 2 . In the region between the curves 
above, 

T = {(x,y)eR 2 >Q | (1/S)x 2 > y > Sx 2 } 



the flow corresponding to the first reaction can go both ways; the sub-tangentiality 

~-2 
1 



condition forces V H T to be parallel to 



In other words, as in Figure 



3.1 



must 



(*) r cannot contain any vertices ojV, and the two sides ofV that intersect 
be parallel to the reaction vector that corresponds to T. 

This is the only requirement that reaction 2X ^ Y (or any other reversible reac- 

tion) imposes on V , since convexity of V guarantees that sub-tangentiality is satisfied 
everywhere else. 

Figure [3T] (b), (c) illustrates the flow corresponding to the two remaining pairs of 
reversible reactions. We may now combine conditions (*) imposed by the three pairs 



of reversible reactions and obtain the polygon V in Figure 3.1 (d) that is forward 



invariant for (3.2). We may realize V from a rectangle in the positive quadrant whose 
sides are parallel to to axes and whose corners are "cut" at intersections with T-type 
regions in the direction of the corresponding reaction vector (see Figure [iP] (d)). The 
original rectangle must be close enough to the axes and "large" enough to contain c(0) 
in its interior and to allow for cuts in regions where the curves y = Sx a , y = (l/S)x' J , 
a € {1,2} are ordered with respect to a. For instance, the cuts at the (0,0) or SW 
corner are in a region where 

8x 2 < (l/S)x 2 <Sx< (l/5)x, 

whereas in the NE corner we have 

Sx < {l/8)x < Sx 2 < {l/5)x 2 . 

This arrangement of the curves avoids the ambiguity of cuts intersecting more than 
one T-type region and is necessary for obtaining a convex V '. 



We therefore conclude that the mass- action system (3.1) is persistent. In fact, 
our informal discussion can be made rigorous and extended to apply to any reversible 
network: 

Theorem 3.1. Any two-species reversible K-variable mass-action system is per- 
sistent. Although our discussion of Theorem |3.1| has been informal, a rigorous proof 



will be given through the more general result of Theorem 6.2 

In the next sections we extend ideas discussed here for the case of endotactic 
networks which we will define soon. For reversible networks, this section revealed 
an important ingredient in constructing V and proving persistence. Namely, when 
c(t) £ V, the aggregate effect of a pair of reversible reactions P ^ P' is to push the 
trajectory towardss the corresponding T region, therefore towardss the interior of V '. 
The "bad" reaction in P ^ P' , i.e. the one that points outside V is cancelled by the 
other, "good reaction". Both reactions are good if c(t) lies on one of the two sides 
parallel to P' — P. Which reaction is good and which is bad comes from the sign of 
binomials of the form 

c p ' - C ■ c p , (3.3) 
where C is some positive constant, in this case C = 5 or C — 1/(5. Analyzing the sign 



of (3.3) for source complexes P and P' is a recurring theme throughout the paper. 



We will refer to this informally as "comparing source monomials up to a constant. 
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4. Endotactic networks. The results we obtained in this paper are applicable 
to endotactic networks, a large class of two-dimensional reaction networks character- 
ized by a simple geometric property. As we will see shortly, the class of endotactic 
networks is larger than the well-known class of weakly reversible networks. We build 
up to the definition of an endotactic network and begin with a few geometric notions. 

Let 7Z = (S,C, 1Z, k) denote a k- variable mass-action system with two species and 
let S = [X, Y) denote its ordered set of species. 

Definition 4.1. By a slight abuse of terminology, the set of lattice points corre- 
sponding to sources of 1Z : 

SCiJZ) = {(m, n) € Z> such that mX + nY G C is a source complex}. 

is also called the set of source complexes oflZ. The source monomial corresponding to 
(m,n) € SC(K) is x m y n . 



Pi 




Fig. 4.1. Essential supports corresponding to vectors vi and V2 

Let A be a set of points and v a vector in the plane. 

Definition 4.2. (i) Let supp v (A)>q denote the minimal half-plane which con- 
tains A and the positive direction of v and is bounded by a line perpendicular to v. 
The v-support of A is the line that bounds supp v (A)>o, 

supp v (A) — dsupp v (A)>q. 

We define supp v (0) = 0. 

(ii) The v-essential subnetwork of a reaction network 1Z is the set of reactions in 1Z 
whose reaction vectors are not orthogonal to v, i.e. 

1Z V = 1Z\{P -+ P' | (P' - P) ■ v = 0} 

(Hi) The v -essential support of a reaction network 1Z is 

esupp v (7Z) — supp v (SC(TZ V )). 

Note that 1Z V = implies that all reaction vectors in 1Z are orthogonal to v. In 
this case we have 1Z V — esupp v (1Z) — 0. 

Example. In Figure |4.1| the dots represent six complexes and the arrows describe 
the four reactions among them. The two lines shown in the picture are the essential 
supports corresponding to the two vectors Vi and V2. 
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Throughout this paper we will also use the notation 

esupp v {TZ) <0 = {P el 2 | (P - Q) ■ v < for all Q £ esupp v (K)} 

and define esupp v (72-) >o similarly. 

Definition 4.3. A reaction network 1Z is called endotactic if for any nonzero 
vector v with 1Z V / w have 

{P^ P' \ Pe esupp v {K) and P' € esupp v (K) <0 } = (4.1) 



In other words, and perhaps more illuminating, a network is endotactic if and 
only if it passes the "parallel sweep test" for any nonzero vector v : sweep the lattice 
plane with a line L orthogonal to v, going in the direction of v, and stop when L 
encounters a source complex corresponding to a reaction which is not parallel to L. 
Now check that no reactions with source on L points towards the swept region. Note 
that if 1Z V — then L never stops. In this case we still say that the network has 
passed the parallel sweep test for v. 

It is easy to verify that a reaction network is endotactic. Although, by definition, 
condition (4.1) needs to be met by any vector v, it is easy to see that only a finite 
number of vectors v need to satisfy (4.1) in order to conclude that the network is 
endotactic. This fact is explained in Proposition 4.1 which we discuss next. 

Let V denote the set of inward normal unit vectors to the sides of conv(SC(7Z)). 
If conv(SC(lZ)) is a line segment, we consider both normal unit vectors as inward. 

Proposition 4.1. A reaction network 1Z is endotactic if and only if condition 
(4-1) holds for any vector v 6 V U {±i,±j} 7 where {i,j} is the standard base of the 
cartesian plane. 

Proof. The "only if implication is clear. For the "if implication, let v be any 
unit vector such that 1Z V ^ 0. If v € V U {±i, ±j}, there is nothing to prove. Else, if 
[Qi ■ • ■ Qi] denotes the polygon dconv(SC(lZ)), assume that esupp v (lZ) = {Qk}- 




Fig. 4.2. Possible reactions Qk — > P' in the proof of PropositionUjJ] In (a), the lines Qk—iQk 
and QkQk+l are essential supports for their inward normal vectors; the target P' belongs to the 
interior of the angle Qk-iQkQk+i ■ In (b), only QkQk+l essential support; P' is constrained on 
the half line [Q^Qk-l- I n both cases, P' £ esupp v (R.) >q . 



First suppose that conv(SC(lZ)) is not contained in a line. Let Vk-i and Vk 
denote inward normal vectors to Qk-iQk and QkQk+i- The source complex Qk may 
belong to both or exactly one of the sets esupp Vk (H) and esupp Vk+1 (TV). 

If Q k € esupp Vk (K) n es-upp Vk+1 (^), then (see Figure [4~2| (a)) condition p"!} 
applied to Vk and Vk+i shows that the target P' of any reaction Qk — > P' must 
belong to the closed positive cone generated by the vectors QkQk-i and QkQk+i- 
This cone (except its vertex Qk) is contained in esupp v (72.) >o- 
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(d) 



Fig. 4.3. Examples of endotactic networks - (a) and (c), and non-endotactic networks - (b) 
and (d). If the convex hull of the sources is a non-degenerate polygon ((a) and (b)), the parallel 
sweep test is done with vectors o/V. If the sources are confined to a line ((c) and (d)), the parallel 
sweep test for {±i, ±j} needs to be performed also; note that using just the vectors in V would have 
concluded that the network (d) is endotactic. 



Vk (lZ)\esupp Vk+1 (TV) , all reactions Qk — > P' are along the line 

ic+i shows 

(b)) p' e e S u P p v (n) >a . 



If Qk S esupp' 

iQk- In fact condition (4.1) applied t o V k+i shows that they must have the 



Qk 

direction of QkQk-v Therefore (see Figure 

Finally, if conv(SC(7Z)) — [Q1Q2} is a line segment then (4.1 ) applied to the two 



4.2 



normal vectors of the line Q1Q2 shows that all target complexes lie on the line Q\Q2- 



Moreover, (4.1) applied to {±i,±j} guarantees that in fact the all target complexes 
of reactions with sources Q\ and Q2 lie on the segment [Q\Q2\- It is now clear that 
(4.1) holds for any vector v. □ 

Using Proposition |4.1[ it is straightforward to check whether a network is endo- 



tactic or not. A few examples are presented in Figure [4~3| 

We have mentioned that the class of endotactic networks includes the class of 
weakly reversible networks. Indeed, this is easy to see: 

Lemma 4.4. Any weakly reversible reaction network is endotactic. 

Proof. In a weakly reversible network all complexes are sources, SC(1Z) = S. 
The definition of the essential support implies immediately that for every vector v all 
reactions originating on esupp v (1Z) point inside esupp 1 (IZ)^. □ 

Because of its relation to complex-balancing, weak reversibility is a widespread 
assumption in much of the work done in Chemical Reaction Network Theory ([S]). 
On the other hand, the results of the present paper naturally lend themselves to the 
much larger class of endotactic networks. In the subsequent sections we will prove 
persistence and permanence of endotactic re-variable mass-action systems with two 
species. We conjecture that the same results are true for any number of species: 

Extended Persistence Conjecture. Any endotactic n-variable mass-action system 
is persistent. 

Extended Permanence Conjecture. Any endotactic n-variable mass-actioi^\ sys- 
tem is permanent. 

Here "extended" refers to two different generalizations: the assumption of an 
endotactic, and not necessarily weakly reversible network, and for the use of the re- 
variable mass-action kinetics, and not the usual mass-action kinetics, where re is a 
fixed vector of parameters. 

5. Construction of a forward-invariant polygon. As in the proof of persis- 



tence of the reversible example (3.1 ), our approach to proving Theorem 6.2 is based 



5 Same for power-law. 
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on defining a convex polygon V whose convex hull contains T(c(0)). Here we explain 
how this construction is carried out for a not necessarily reversible network. 



In example (3.1) the aggregate effect of any two reverse reactions P ^ P' is to 
drive the trajectory towards the T region 5 < c p ~ p < 1/5, therefore towards the 
interior of V . In other words, the source of the "bad" reaction that might force the 
trajectory to cross to the exterior of V is smaller (up to some constant) than the 
source of the "good" reaction that keeps the trajectory in the interior of V . 



The proof of persistence for example (3.1) contains all the ingredients for the 
proof of the persistence for any two-species reversible K-variable mass-action system. 
In the more general setup of a weakly reversible network, we need to adjust the setup 
of the invariant polygon and parts of our discussion. 

In the reversible case we saw that if c(to) € V, every pair of reactions P ^ P' 
contains a "good" reaction that dictates the overall effect of the reversible reaction. 
In the weakly reversible case we show that there is a "good" reaction that dictates 
the direction of c(to)j i-e. the aggregate effect of the whole network, up to a small 
enough perturbation. This amounts to showing that at time to with c(fo) € V there 
is a complex that is source to a "good" reaction and whose corresponding monomial 
is larger than the monomials corresponding to all other source complexes (up to some 
constant). This way we are led to comparing monomials c p and c p for all pairwise 
combinations of source complexes P and P'. 



Let us fix a trajectory T(c(0)) of (2.6) and denote c(t) = (x(t),y(t)) for t > 0. 



Let (1, r\) . . . , (1, r e ) and (1, s\) . . . , (1, s/) be normal vectors of edges with finite and 
nonzero slope in the complete graph with vertices SC(1Z) : 

{ri, . ,.,r e ,si, . .. ,Sf} = I — — — | (mi,ni), {m 2 ,n 2 ) € SC{TZ),m\ ^ m 2 ,n 1 ^ n 2 

I n 2 ~ni 

where s± < . . . < s / < < r\ < . ■ . < r e . 

Let < 5 < 1 be fixed. The pairwise comparison of source complexes P = (m, ri) 
and (m! ,n') up to a constant leads to c p — Sc p — and c p — (l/<5)c p = 0, or 
y = tfVC"-™')^ and y = (1/ l S) 1 ^ n ~ n '^x' T , where a G {n, . . . , r e , s u . . . , s/}. Let 

5' = mm{S 1/{n '- n) | (m,n),(m',n') G SC(K),n ri) 
and consider the curves 

y = 5'x r * ,y = (l/S')x Ti and y = S'x Sj ,y = (\/5')x Sj for all i G {1, . . . ,e} and j G {1, . 



depicted using dots in Figure 5.1 We further let r i+ 5 be positive reals and Sj+.g be 
negative reals for all i G {0, . . . , e} and j G {0, . . . , /} such that 

ri_. 5 < ri < r i+ , 5 and Si_. 5 < s l < s i+ , 5 . 

The solid curves in Figure |5.1| are given by 

y = x Ti + - 5 and y = x Sj+ 5 , i G {0, . . . , e} and j G {0, . . . , /}. 

Remark 5.1. In our construction, the vertices of V will be confined to these 
"fractional indices" curves. Although this constraint is not necessary and we may 



allow more freedom in choosing the vertices of V as in the reversible example (3.1), 
our construction conveniently makes V depend continuously on a single parameter. 
Let £ < 1 be positive and small enough and M > 1 be positive and large enough 



such that (see Figure 5.1 1 
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Fig. 5.1. Setup of the invariant polygon V = A\ . . . A e +iBi . . . Btj r \G\ . . . C e +\D\ . . . -D/+i 



2. 



(PI) c(0) € (£, M) 

(P2) all intersection of the curves in Figure 5.1 lie in the square (£,M) 2 ; 
(P3) (0,£) 2 and (M, oo) 2 lie below, respectively above all negative-exponent 
(sj or sj + , 5 ) curves. Likewise, (0,£) x (M, oo) and (M, oo) x (0,£) lie above, 
respectively below all curves of positive exponent (r^ or r^+.s); 
(P4) All negative exponent curves intersect the line segments (0,£) x {M} and 

{M}x(o,0; 

(P5) £ < min(<5 1/( ™'- n) | (to, n), (to', n') € SC{H),n ^ n'), 

( < min( ( 5 1/(m '- m) | (m, n), (to', ra') € SC(ft), m ^ to'), 
M > max((l/^) 1/( "'- n) | ( m, n), (to', n') € SC(1Z), n ^ n'), 
M > max((l/ ( 5) 1 /( m '- m ) | (to, n), (to', n') e SC(TZ), to ^ to'). 

The goal is to obtain V = Ai . . . A e+ iBi . . . Bj +1 Ci . . . C e+ \Di . . . Df + i as in 
Figure |5.1| with the property that 

(P*) the dotted curves y — S'x' 7 and y = (l/(5')a; CT intersect the two sides of V that 
are orthogonal to the vector (1, a) and only these two sides. 

Conditions (P1-P5) allow for the construction of such a polygon, and moreover, 
assure that V is invariant, as we will see in the proof of Theorem 6.2 Due to (P4), the 
Bi and Di vertices are easily constructed in the correct shaded regions. Condition (P2) 
guarantees that on the shaded areas the pair of dotted curves y = 8'x a , y = (l/<5')x CT , 
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(CO) 




(0,0 



Fig. 5.2. Detail of vertices Ai, 



a £ {ri, . . . , r e , Si, . . . , s/} lie between the corresponding solid curves, an important 
fact that we will take advantage of later. Conditions (P3) and (P5) will also come in 
handy, as it will make some pairwise comparisons of source monomials straightforward. 
Finally, condition (PI) make the initial condition lie in conv(P). 



We construct V = A 1 . . . A e+l B l . . . B f+1 Ci . . . C e+1 D 



D 



one vertex at 



the time starting with Ai = (a,a s - 5 ) and going counterclockwise. All vertices of V 
lie on a fractional index curve and satisfy property (P*) above. In other words the 
sides AiA i+ i, dCi + i are orthogonal to (l,rj) and the sides BjBj +1 and DjDj + i are 
orthogonal to (l,Sj). We also require that A e +\B\, C e+ \D\ are horizontal and that 
Bf + \C\ is vertical. The polygon V constructed this way is unique. 

The vertices of V can indeed be chosen to lie in the four shaded areas of Figure 



5.1 



For instance, lim Q _i.o Ai — (0, 0); therefore Ai £ (0, £) 2 for all i € {1, . . . , e + 1} 
for small a, whereas (P4) guarantees that Bj £ (M,oo) x (0, £) for all j. Finally, note 
that the side Df+iAi of V need not be vertical, but in the limit a 0, Dt+iAi is 
vertical. In particular, Via) is convex for small a. 

Let ao be such that V(oiq) is a convex polygon whose construction is as explained 
above. Then a € (0, ao] defines a one-parameter family of convex polygons; the 
generic V{a) is outlined in Figure 5.1 The vertices A\, . . . , A e are illustrated in more 



detail in Figure 5.2 



Finally, let us remark that P(a) varies continuously with a and conv(V(a)) is 
decreasing in a € (0, ao) and is a cover of M.y Q : 



conv(V(a')) D conv{V(a")) if a' < a" and conv(V(a)) = 



a 2 



(5.1) 



6. Endotactic two-species K-variable mass-action systems: persistence 
and permanence. 

6.1. Persistence. Let (TZ, S,C,k) denote a two-species, endotactic K-variable 
mass action system with Kp^p>(t) £ (r), l/f]) for all P P' € TZ. Recall that c(t) = 
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(x(t),y(t)) denotes a solution of (2.1) with initial condition c(0). 

The following lemma is a key ingredient in the proof of Theorem |6.2| It states 
that if for some to, the source monomial of the reaction Po — > Pq is larger than all 
other source monomials (up to a constant), then c(to) does not deviate too far from 
Po" -Po- 

Lemma 6.1. Let P — > Pq G TZ and let v be a vector such that (Pq — P ) • v > 0. 
Also let U C SC(lZ)\{Po}. There exists a constant 5 such that if for some t > we 
have c(to) € K?, and 

c{to) P < Sc{to) Po for allP eU 

then 

' Kp ^ PL {to){P'o ~ Po)c(t ) P ° + Yl *P^P>(to)(P' - P)c{tof ] ■ V > 0. 

Proof. We take 



8 = 



^(P^-Po)-v 



Ml £ \\p'-p\\ 

and we have 

K Pa ^P>(to)(Po - PoMtof + *P^P'(to)(P' - P)c(to) P I • v > 



p^p\Peu 



>Kp ^(to)c(to) Po (Po-Po)-v~ J2 («P-i"(*o) l|v||||P'-P|| Sc(t ) Po J > 

p^P',Peu \ ) 

>UPo-Po)-v-5(\\v\\/ v ) I|p'-p||Wo) Po = o 

\ P^P'GK ) 

□ 

We may now state and prove our main result. 

Theorem 6.2. Any two-species, endotactic n-variable mass-action system is 
persistent and has bounded trajectories. 

Proof. Recall that V denotes the set of inward normal vectors to the sides of 
conv(SC(lZ)). Let n e V U {±i, ±j}, and let P n — > P^ denote a fixed reaction such as 
P n € esupp n (lZ) and P^ £ esupp n (lZ) >0 . Since 7?. is endotactic, such a reaction exists 
for all n eVU{±i,±j}. 

Note that including {±i,±j} here guarantees that at least one vector n £ V U 
{±i,±j} satisfies esupp n (lZ) ^ 0. For every such vector n we let 5 n denote the 
constant from Lemma |6.1| that corresponds to reaction P n — > P„, v = n and U = 
SC(K) n esupp n (K) >0 and 

S = min <5 n . (6-1) 

n 

Then 6 < 1 and the construction from section [5] yields the convex polygon 

V(ao) = A\ . . . A e+ iP>i . . . Bf+xCi . . . C e +iDi . . . P/+i = V for simplicity of notation. 
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We view V as the union of its sides, so that V — d(conv(P)). By construction, c (0) g 
(£, M) 2 c conv{V). We show that conv(V) is forward invariant for the dynamics (]2.6[), 



i.e. c(t) £ conv{V) for all t > 0. More precisely, we prove that c(i) cannot cross to 
the exterior of V by showing that if c(to) £ V for some to > 0, then 

c(i ) • n > whenever c(t ) e 7> (6.2) 

where — n is a generator of the normal cone N-p(c(to)) of V at c(to). Recall that 

N v (c{t )) = {v £ R 2 | v • (x - c(t ))) < for all x £ conv(P)}. 

Note that n is orthogonal to one of the sides of V and points inside V . 



Inequality (6.2 1 is rewritten as 



v p->P'eK, 



Kp ^ P ,(t ){P' -P)c(toY 



n > 



(6.3) 



and is satisfied if lZ n = 0. Otherwise, let P n —> P^ £ lZ n with P n £ esupp n (lZ). 



We rewrite the left hand side of (6.3) by separating the reactions with source on the 
n - essential support of 1Z and further isolating the effect of the reaction P n -» P^ : 



( J2 Kp^P>(to)(P'-P)c(t ) P )-n = 

\P^P'£R„ / 
= ( E KP^P>(t )(P' -P)c(t ) P + 



P£esupp n (TZ) 

+Kp^p k (t )(P^- Pn)c(t ) P " + K P ^ P >(t )(P' -P)c(t ) P \ -n. 

P^P'e1Z n ,P$:esupp"-(1l) } 

Since all source complexes of lZ n lie in esupp n (lZ I1 )>o, the reaction vector P' — P with 
source P £ esupp n (7Z) satisfies (P' — P) ■ n > 0. Note that for the reaction P n — > P^ 
the inequality is strict, 

(i* - P n ) ■ n > Q. 
It is therefore enough to show that 



KP»-^(*0)(i^ - ^n)c(t ) P " + Yl KP^P'(t )(P' - P)c(toY 



V 



P(£esupp n (7Z) 



J 



n > 



(6-4) 

in order to verify (6.3). In turn, (6.4) will follow from Lemma 6.1 with li = SC(1Z) D 



esupp n (lZ) > o and the fact that 

6c(t ) p " > c(t f (6.5) 
for all P £ U, inequality whose verification will complete the proof of the theorem. 
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T77I ~ n 




Fig. 6.1. Positions of a source monomial P relative to P z 



The proof of inequality (6.5) depends on the position of c(to) on V . More pre- 
cisely, c(to) might belong to one of the four polygonal lines [Ai . . . A e+ i], [B\ . . . 
[Ci . . . C e +i] and [D\ . . . -D/+i] or to one of the four remaining sides of V . We discuss 
the first case, c(i ) G [Ai...A e+1 ]. Suppose c(i ) = {x(t ), y(t )) lies on the side 



[j4iAi + i]. Then (see Figure 5.2) we have 



(l/<5'MioP +1 <y(to) <6'x{to) r '-\ 

where only one side of the inequality holds if i = 1 or i = e. 

Let P = (m,n) and recall the notation P n = (m n ,n n ). We let a = ™~^* if 
n n 7^ n. Relative to P n , the source complex P may belong to one of the three shaded 
regions in Figure |6.1| 

Region I. Here m > m n and n < n n . Also, a > ri and therefore a > r^+i and we have 

y(to) > (l/<5')^(^o) r ' +1 > {l/Sfl^^x^y, 
which implies 8x{to) m ^y{t Q ) n ^ > x(t Q ) m y(t ) n . 

Region II. If m ^ m n and n ^ n n we have m > m n , n > n n and cr < 0. Condition 
(P3) implies 

i/(t ) < s'x(t y < s'/^-^xitor, 

which again implies x(to) m y(to) n < i5a;(to) m "y(^o)" n - If ™ = TO n or n = n n we need 
to show 6 > or 5 > x m ~ mn , which are immediate from (P5). 

Region III. Here m < m n and n > n n . We have a < Vi, so cr < rj_i and 

yfa,) < S'x(t p-i < 5 l l^- n ^x{t y, 

therefore 5x(t ) m "y{t ) n " > x(t ) m y(t Q ) n . 

The cases when c(t ) lies on the other three polygonal lines [B\ . . . \C\ . . . C e+ \] 

and [D\ . . . -D/+i] arc similar. Finally, suppose that c(to) is on one of the four remain- 
ing sides, [A e+ iSi], [B/+iCi], [C e+ i£>i] and [D/ + iAi], for instance c(t ) e [A e+ ii?i]. 
Then we must have n > n n and 

y(t ) < S'x(t y < S'/^-^xitor, 

or x(t ) m y(t ) n < 6x(t Q ) m "y(t Q ) n ". □ 

Remark 6.1. Note that if the inclusion n(t) € (q^l/rff 1 is not necessarily true 



for all t > 0, but holds for t = to then the inequality (6.2) is true (at time to) 



Kp^P'(t )(P' -P)c(t a ) p n>0 (6.6) 



\P^P l en 
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Lower-endotactic networks. A useful observation can be made from the proof 
of Theorem |6.2| Namely, we need not require that the network is endotactic in order 
to conclude that a bounded trajectory T(co) is persistent. Instead, it is enough that 
the network is merely lower-endotactic, meaning that it passes the parallel sweep 
test for vectors that point inside the closed positive quadrant. Let be the set of 
inward normal vectors to the sides of strictly negative slope of SC(1Z) (see Proposition 



4.1 1. In other words, a network is lower-endotactic if condition ( |4.1[ ) is satisfied for 
all vectors v € U {i,j}. Indeed, we may construct the invariant polygon V such 
that the square [0, M] 2 contains T(c ). Then, if c{t) reaches V, it must be on the 



polygonal line [D f + \A\A2 . . . A e+ iB{\. and the proof of (6.2) only requires that H. is 
lower endotactic. We cast this result as a corollary to the proof of Theorem |6.2| 

Corollary 6.3. Any bounded trajectory of a lower-endotactic K-variable mass- 
action system with two species is persistent. 

6.2. Permanence. The proof of Theorem |6.2| allows us to say more about the 



dynamics of (2.1). Indeed, without much additional effort we can prove the following 



stronger version of Theorem 6.2 



Theorem 6.4. Any two-species endotactic K-variable reaction system is perma- 
nent. 

Proof. Suppose that for some small e > and for some a > we can construct 
the polygon ^(ao + e). Recall that the polygon V = V(a) depends continuously on a 
and the function conv(V) : (0,ao + e) — !• IR is strictly decreasing, i.e. conv(V(a)) D 



conv(V(a')) whenever a < a' < ao + e. Let T(c(0)) denote a trajectory of (2.1 ) with 
initial condition c(0). We will show that the trajectory T(c(0)) eventually ends up in 
conv(V(a )). 

Since U Q e(oa +e] conv{V(a)) = R^q we may choose < ct\ < ao such that 
0+e ] conv(P(a)). Let $ : \J a& [ aitOl0+e ]'P(a) — ► [ai,ao + e] be defined 
as <&(x, y) = a if (x,y) € V(a). We need to show that $(c(i)) > a for large enough 
t. Suppose this was false; then, since < I>~ 1 [q!o, oo) = conv(V(ao)) and $ _1 [ai,oo) = 
conv(V(ai)) are both forward invariant sets for (2.1), we must have that <I>(c(£)) € 
[ai, ao] for all t > 0. 

Note that <& is differentiable on its domain except at the points on the fractional 
index curves y = x a for all a G {r.5, . . . , r e _|_.5, s.5, . . . ,Sf + 5 }. Let c € IR^q be such a 
point, and let $1 and $2 denote the smooth functions that define $ on the two sides 
of the fractional curve in a neighborhood of c. The subgradient [501 Definition 8.3] of 
$ at c is 



S$(c) = (aVfi(c) + (1 - a)V$ 2 (c) | a £ [0, 1]} (6.7) 

and therefore <I> is strictly continuous [20, Definition 9.1]. Note that ( |6.7[ ) is also 
valid for points c that do not belong to a fractional index curve and in that case 
V$i(c) = V$ 2 (c) = V$(c). Since c(t) is smooth, it follows that $ o c(t) is also 
strictly continuous; in particular, a generalized mean value theorem [201 Theorem 
10.48] implies that for all i > there is r t G [0,t] such that 

<E>(c(i)) - $(0(0)) = s t t for some scalar s t € <9($ o c)(r t ) (6.8) 

The chain rule for subgradients [201 Theorem 10.6] 



<9($ o c)(t) C {v • c(t) I v e 0$(c(t))} 



(6.9) 
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connects the subgradient of $ o c to that of $. Note that the proof of Lemma [6.1| 
can be modified to show that for any compact set K C R> there exists £ > such 
thatV$i(e(£)) • c(i) > C and V$ 2 (c(t)) ■ c(i) > C for c(t) £ X. This, together with the 



special form (6.7 1 of and together with (|6.9| shows that 



inf 5($oc)(t) > C 



From (6.8 1 we obtain 

$(c(t)) > $(0) + Ct for all t > 0, 
which contradicts the fact that <£>(c(i)) £ ao] for all i > 0. □ 

Remark (existence of positive equilibria for two-species autonomous endotatic mass- 
action systems). The existence of a positive equilibrium has been shown for weakly 
reversible mass-action systems in [7] . For endotactic networks we may observe that if 
the rate constant function depends on both time and the phase point c, k = n(t, c) £ 
(?/, 1/?/) then the set of all possible trajectories is still the same as for n — n(t), and 
our persistence and permanence results still hold. In particular, if n — n{c) £ (77, 1 / 77) 

together with the Brouwer Fixed Point Theorem 



6.4 



for all c £ K>o> then Theorem 

guarantees the existence of a positive equilibrium in conv(V(ao)). 

7. The Global Attractor Conjecture for three-species networks. As 

mentioned in the Introduction, the Global Attractor Conjecture concerns the asymp- 
totic behavior of complex-balances mass-action systems and is of central importance 
in Chemical Reaction Network Theory. 

Definition 7.1. A point c* £ ]R™ of a mass-action system (1Z,S,C, k) is called 
complex-balanced equilibrium if the aggregate flow at each complex oflZ is zero. More 
precisely, for each Pq £ C we have 

P-vPo Po^P 

A complex-balanced system is a mass-action system that admits a strictly positive 
complex-balanced equilibrium. 

A complex-balanced system is necessarily weakly reversible [5] . A very important 
class of complex-balanced mass-action systems is that of weakly reversible, deficiency 
zero systems 019]. The deficiency of a network is given by n — I — s, where n, I and s 
denote the number of complexes, the number of linkage classes and the dimension of 
the stoichiometric subspace. Computing the deficiency of a network is straightforward, 
as opposed to checking the existence of a complex-balanced equilibrium. 

It is known that a complex-balanced system admits a unique positive equilib- 
rium cx; in each stoichiometric compatibility class £ and this equilibrium is complex- 
balanced [8]. Moreover, each such equilibrium admits a strict Lyapunov function 
which guarantees that cs is locally asymptotically stable with respect to E [Sj H2] . 
Recall from the Introduction that the Global Attractor Conjecture states that Cs is 
in fact globally asymptotically stable with respect to E. 

It is known that for complex-balanced mass-action systems, all trajectories with 
positive initial condition converge to the set of equilibria |18j . As a consequence, 
showing that the unique positive equilibrium in the stoichiometric compatibility 
class E is globally attractive amounts to showing that no trajectories with positive 
initial conditions have w-limit points on the boundary of E. 
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In particular, any condition that guarantees the non-existence of boundary equi- 
libria for a stoichiometric compatibility class £ implies that is globally attractive. 
For instance, it was shown that a face F of £ contains boundary equilibria only if the 
species corresponding to coordinates that vanish on F form a semilocking set ( siphon 
in Petri nets terminology) [IJ [2]. It follows that if no face of £ corresponds to a 
semilocking set then £ does not have boundary equilibria and ce is globally attractive. 
Algebraic methods for computing siphons have been devised in |17| . 

Some progress has been made for the remaining case where boundary equilibria 
cannot be ruled out. It was shown that w-limit points cannot be vertices of £ 0311] 
and that codimension-one faces of £ are "repelling" [2]. As a consequence, the Gobal 
Attractor Conjecture holds for systems with two-dimensional stoichiometric subspace 
[2]. For systems with three-dimensional and higher-dimensional stoichiometric sub- 
space the conjecture is open. In particular, previous results do not apply to the general 
three-species systems. In what follows we show the conjecture to be true in this case. 

Theorem 7.2. The Global Attractor Conjecture holds for three-species networks. 

Proof. Let (S,C,7Z, k) denote a complex-balanced reaction system with ordered 
set of species {A, Y, Z} and denote the corresponding concentration vector at time t 
by c(t) = (x(t),y(t), z(t)). Let c € R> be arbitrary and let T(c ) denote the forward 
trajectory of the concentration vector with initial condition cq. We want to prove that 
T(co) converges to the unique positive equilibrium in the stoichiometric compatibility 
class of c . 

The existence of the Lyapunov function guarantees, on one hand, that there is a 
neighborhood of the origin that is not visited by T(cq), and on the other hand, that 
T(co) is bounded [5JQ]. In other words, there exists e > such that x(t)+y(t) + z(t) > 
3e and x(t), y(t) and z(t) are all smaller than 1/e at all times t > 0. 

Our previous discussion explained that it is enough to rule out boundary w-limit 
points for T{cq). We will construct a compact set K. C R> such that T(cq) C JC. 

If TT xy denotes the projection onto the coordinate coordinate plane xy, then 



x(t) 

m 



^ K P ^ P ,z{t) p * (x(t), y(t)r^ p \ xy (P' - P). (7.1) 



where Pz denotes the stoichiometric coefficient of the species Z in P and ir xy (P) is 
the complex with species {A, Y} obtained by removing species Z from P. 

Let Kmin = min{K, 1/k \ k is a rate constant of some reaction in 1Z} and let s max 
denote the maximum stoichiometric coefficient of any species in a complex of 1Z. Let 

f] = K min£ mam (7.2) 

and note that r\ < 1. 

The reactions of 1Z project onto reactions among complexes n xy (C) with species 



X and Y. Let n xv (1Z) denote the resulting reaction network. Motivated by (7.1) we 
consider the k- variable mass-action system of variables x and y given by 



x(t) 

m 



£ ^ P ,{t)(x(t). /y {t)) p {P' -P). (7.3) 



where n{t) 6 (77,1/77) for all t > 0. Since 1Z is weakly reversible, ir xy (7Z) is weakly 
reversible as well. In particular, n xy (TZ) is endotactic and we may construct an in- 
variant polygon V xy for (7.3| as in the proof of Theorem 6.2 Section [5] contains the 



details of the construction and we use its notations in what follows. 
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(d,l/e,l/e) 



Fig. 7.1. The subset of dK that can be reached by T(co). The phase point c(to) may belong to 
one of the three flat parts where one of the coordinate is d, or on one of the cylindrical parts along 
the axes. 



In Section [5] the vertices of V lie on the "fractional indices" curves to make 



V = V(a) vary continuously with a. As noted in Remark 5.1 this constraint is not 
necessary for proving that V is an invariant set for the corresponding K-variable mass- 
action kinetics. For the purpose of this section we will instead construct V = V xy 
such that A\Df + \ is a vertical line and that its distance to the y axis is equal to the 
distance from the horizontal line A e+ iBi to the x axis; let d > denote this distance. 
Since V xy may be chosen to contain any compact subset of R 3 (see (5.1 )), we assume 
that it satisfies (using the notations from Figure 5.1): 



(7.4) 



We construct in the same way V yz and V z 



the square [e, 1/e] 2 is included in the square [£, M] 2 . 

and we define 

JC = {(x,y,z) € [0, l/e] 3 | (x,y) £ conv(V xy ), (y, z) € couv(V yz ), {z,x) G conv(V zx )}. 

(7.5) 

Note that JC is a compact subset of IR^g and Co G JC. We show that the trajectory 
T(co) is included in JC. To this end, suppose that at some time point to the trajectory 
reaches the boundary of JC. It is then enough to check that 



n • c(t ) > 



(7.6) 



where — n 6 N/c(c(to)) is any generator of the normal cone of JC at c(to). 

Let C xy de note the subset of V xyi given by the polygonal line D f +1 AiA 2 . . . A e+1 Bi 



(see figure 5.1 1 and similarly define C yz and L zx . Because c(t) € (0, 1/e) 3 for all t > 0, 



(7.4) implies 

c(t ) £ [0, 1/r?] 3 n ((£ xy x K >0 ) U (C yz x M >0 ) U {C zx x M >0 )) (7.7) 
Since x(to) + y(to) + z{to) > 3e, at least one of x(to), y(to) and z(to) is larger than 



e. Assume z(t ) > e. Then (7.7) implies that 



c(t ) e v, 



■>■!) 



(7.8) 



and therefore the z coordinate of n is zero. It follows that inequality ( 7.6 1 is equivalent 
to 



7Txy(n) 



Kp^P>Kto) P Hx(t ),y(t )y™ (P \n X y(P')-n xy {P))}>0. (7.9) 
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where Pz denotes the stoichiometric coefficient of the species Z in P. Note that 
n X y (n) is a generator of the normal cone of V xv at (x(t ),y(t )). Our c hoic e of 77 



( 7.2 1 implies that np^pi z(to) Pz £ (77,1/77) and (7.9) is implied by Remark 6.1 □ 



8. Examples. 

8.1. Thomas-type models. The Thomas mechanism (see [T31 ch. 6]) is a 
substrate inhibition model for a specific reaction involving oxygen and uric acid in 
the presence of the enzyme uricase. After nondimensionalization the ODEs for oxygen 
(v) and uric acid (u) become 



du 
~dl 
dv 
~dt 



= a — u — T(u, v)uv, 
= a(b — v) — T(u, v)uv, 



(8.1) 



where in the Thomas model T{u, v) 
and K are positive. Using Theorem 6.4 we can show that 



P(l 



Ku 2 ) 1 . Here all parameters a, b, a, 



.1) is permanent for 



any function T(u, v) > that is continuous and does not vanish on a compact K D 
[0, a] x [0,6] (in particular, this is true for function T in the Thomas mechanism). 
Indeed, for some time point to we have (u(t),v(t)) € K] continuity of T guarantees 
that for t > to we have 77 < T(u(t),v(t)) < 1/rj for some 77 > 0. The dynamical 
system (8.1) can be written as a k- variable mass-action system with reactions given 
(a), where the reaction rates are specified on the reaction arrows and 



.1 



in Figure 

T(t) = T(u(t),v(t)). This network is endotactic and Theorem |6.4| implies that the 



dynamical system (8.1) is permanent 



8.2. Power- law systems. In the proof o f Th eorem 6.2 we have not used the 
fact that the exponents P of monomials c p in (2.1) arc nonnegative integers; this is 



the case when P C Z?, represents a chemical complex, which has been our framework 



thus far. In fact, the proof of Theorem 6.2 accommodates the case when P C K 2 has 
any real components. 

More precisely, keeping the notation c(t) = (x(t),y(t)), we consider an ODE 
system of the following form: 



Pesc vev(P) 



(8.2) 



where SC C M 2 is the set of "sources" P and V(P) is the set of "reaction vectors" 
with source P. We also suppose that Kp tV (t) G (77, 1/rj) for some < 77 < 1. With this 
interpretation, (|8.2[ ) is a generalization of a K- variable mass-action ODE system (2.1 ) 
and Theorems |6.2| and |6.4| still apply. 



A particularly important example of power-law, not necessarily polynomial sys- 
tems is the class of S-systems (15 ), where each component of the right hand side 
of ( |8.2| consists of a difference of two "generalized monomials" (i.e., monomials with 
real exponents). S-systems are common in the modeling of metabolic and genetic 
networks. For example, consider the following S-system: 



dx 

~dt 
dy 

d.t 



2x- 1 y 1 - 5 -y 
y~ 2 — Vbx~ 1 y 



(8.3) 
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(a) u 





Fig. 8.1. Examples of reaction networks: (a)-Thomas model \8. 1\ ; (b)-the S-system 
(c)-Lotka-Volterra system \8.5\. 



Note that it is not obvious that trajectories of (8.3) cannot reach 9IR> in finite 



time. However, using Theorem 6.4 we can easily see that (8.3) is in fact permanent 



and, in particular, all trajectories with positive initial condition will be bounded 
away from zero and infinity. Indeed, the generalized monomials in (8.3), i.e. the 



points (—1, 1.5), (0,0.8) and (0, —2), as well as the co rres ponding "reaction vectors" 



(2, — y/5), (— 1, 0) and (0,1) are illustrated in Figure 
endotactic and Theorem |6.4| applies. 



.1 



(b). This configuration is 



Note that the same conclusion holds for any system of the form 



dx 
~dt 
dy 
dt 



2K 1 (t)x~ 1 y 1 - 5 



K 2 (t)y 



0.8 



(8.4) 



where Ki(t) G (r], l/rf) for all i € {1, 2, 3} and for all t > 0. 

8.3. Lotka-Volterra systems. The classical two-species predator-prey model 



A^2A A + B^2B B^>0 



(8.5) 



is neither endotactic nor lower endotactic. Since, for fixed parameters its trajectories 
are either constant or closed orbits, the system is not permanent (see Figure 8.1 
(c)). On the other hand, the fixed-parameter system has bounded trajectories and is 
persistent, which seems to suggest that the requirement of entotactic network may be 
weakened in Corollary |6.3[ Note, however, that the result of the corollary concerns 
K-variable mass-action systems. It is not hard to show that, in general, the n-variable 



Lotka-Volterra system (8.5) is not persistent 



8.4. Examples for the three-species Global Attractor Conjecture. We 

present two examples of reaction networks for which no previously known results can 
resolve the question of global asymptotic stability, but for which Theorem 1 7 . 2 1 applies . 



(a) A^B^A+B^A+C 



(6) A + B^A + C A- 



->B 



2A 



(8.6) 



The first example, given by ( |8.6[ ) (a) is Example 5.4 from [2]; it is reversible with 
one linkage class. On the other hand, example ( |8.6[ ) (b) is weakly-reversible, but not 
reversible and has two linkage classes. Moreover, the deficiency of both networks is 
zero and therefore the two mass-action systems are complex-balanced. 
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In both examples all stoichiometric compatibility classes are equal to M> . If 
(x(t),y(t), z(t)) denotes the concentrations of (A, B,C), then it is easy to see that all 
points on the nonnegative z axis are equilibria for both examples. Therefore boundary 
equilibria exist, and moreover, except for the origin, they all lie on a codimension-two 
face of R>q. In this case all previously known results stay silent, but Theorem 



7.2 



guarantees that, in each example, the unique positive equilibrium is globally asymp- 
totically stable. 
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